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^ ! ABSTRACT 
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Context. We present near-IR spectroscopy of a sample of luminous (Mb < -27.5; Lbd > lO'"* Lq), sub-millimeter-detected, dusty (Mj ~ 
10' Mq), radio-quiet quasi-stellar objects (QSOs) at z ~ 2. 
^ Aims. A primary aim is to provide a more accurate QSO redshift determination in order to trace kinematics and inflows/outflows in these 
sub-mm bright QSOs. Additionally, the Ho- and continuum properties allow an estimation of the black hole mass and accretion rate, offering 
insights into the starburst-AGN connection in sub-mm bright QSOs. 
^ ' Methods. We measure the redshift, width, and luminosity of the Ha line, and the continuum luminosity near Ho-. Relative velocity differences 
. between Ha and rest-frame UV emission lines are used to study the presence and strength of outflows/inflows. Luminosities and line widths 
' are used to estimate the black hole masses, bolometric luminosities, Eddington fractions, and accretion rates; these are compared to the star- 
formation-rate (SFR), estimated from the sub-mm derived far-infrared (FIR) luminosity. Finally our sub-mm-bright QSO sample is compared 
] : with other QSO samples at similar redshifts. 
. . ' Results. The Ha emission line was strongly detected in all sources. Two components - a very broad (>5000 km s"') Gaussian and an 
^ ' intermediate-width (>1500 km s"') Gaussian, were required to fit the Ha profile of all observed QSOs. Narrow (<1000 km s"') lines were 
not detected in the sample QSOs. The rest-frame UV emission lines in these sub-mm bright QSOs show larger than average blue-shifted ve- 
rS locities, potentially tracing strong - up to 3000 km s"' - outflows in the Broad Line Region. With the exception of the one QSO which shows 
exceptionally broad Ha lines, the black hole masses of the QSO sample are in the range log Mbh = 9.0-9.7 and the Eddington fractions are 
" " " between 0.5 and ~ 1 . In black hole mass and accretion rate, this sub-mm bright QSO sample is indistinguishable from the Shemmer et al. (2004) 
optically-bright QSO sample at z ~ 2; the latter is likely dominated by sub-mm dim QSOs. Previous authors have demonstrated a correlation, 
over six orders of magnitude, between SFR and accretion rate in active galaxies: the sub-mm bright QSOs lie at the upper extremes of both 
quantities and their SFR is an order of magnitude higher than that predicted from the correlation. 
Conclusions. 

Key words, high-redshift; quasars: emission lines; quasars: general: sub-millimeter: galaxies ; galaxies: active 



1. Introduction 

There is now general agreement that dormant black holes, be- 
lieved to be QSO 'relics', a re common, perhaps ubiquitous, in 
nearby galactic nuclei (e.g.. lSoltanLll982tlFerrarese & Merritt . 



Send offprint requests to: G. Orellana 



2000 ; Tremaine et al. 


20021 Yu & Tre 


mainel 


2002 


Marconi et al. 


, 2004 Merioni & Heinzl 20081 


Shankar , 


2009 



and references therein). The correlation between the mass 
of the black hole and the structural p arameters of the galaxy 



spher oid (Mj,/; ^ 6 x 10 ^ M^nh, e .g. 



19951: iMagorrian et al. 



1998 



Kormendv & Richstone , 



Gebhardt et al 



2000 



2 



Gustavo Orellana et al.: Sub-millimeter Detected z ~2 Radio-quiet QSOs 



Marconi & Huntl 2003 


; Harine & Rix , 
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2009b sueeests 



an intimate connection between the existence of a luminous 
QSO phase and the formation of massive, early-ty pe g alaxies . 
The peak in the space de nsity of QSOs at z~2 d Schmidt et all 
>91t IWall et all l2005h suggests that these massive black 
holes formed at early epochs. It seems that all sufficiently 
massive bulges undergo a QSO phase with the AGN fueled 
by the initial collapse of the host object at high redshift. This 
picture is suppoited by the detection of massi ve elliptical 



galax i es hosting all luminous, nearby QSOs dBahcall et al. 



19971: lMcLureetallll999h 



Any account of bulge formation 
disregards the coevolution of the central 
vice versa- is fundamentally incomplete 



evolution which 
black hole - or 

jKauffmann & Haehnelti I2OO0I: IShankaii l2009b . A continuing 
task is to unify different galaxy types with varied luminosities 
into a common evolutionary scenario or parameter space. 

Over the last few years we have pursued a highly 
successful program aimed at elucidating the star forma- 
tion history of op t ically-bright radio-q u iet QSO host galax- 
ies ( Pridd ev et al, ,_ lDr.i-.cr.n ot ai ' 
2003al: lO mont et al 
2002; Omont et al.. 



20071 iR obson et al.'. '2004^; 'Priddev et al.', 
IT 



2003t Priddev et al. 



2003b; Isaaketal., 
2001; McM ahon et al.L 1999; O mont et all 
1996a.: .McMahon et ali .1994 llsaak et al.L (1994) via the de- 
tection of thermal emission from cool (~50 K) dust with 
SCUBA, MAMBO and their single-pixel/few-pixel predeces- 
sors. Optically-bright QSOs were targeted as likely signposts 
to massive underlying host galaxies in an era before large blind 
sub-mm surveys were possible. This program concentrated 
on radio-quiet (to minimize the contamination of synchrotron 
emission to the sub-mm fluxes) QSOs at both the highest 
redshift (z>4) and at the epoch of peak QSO activity (z~2). 
Typical star formation rates (SFRs) - estimated from the con- 
version of sub-mm fluxes into far-infrared (FIR) luminosities 
and then SFRs (e.g., 



McMahon etallll999h - of ~ 500 MJyv 



were found in sub-mm detected QSOs, although rates upto 
~ 2500 Mo/yr are observed in some cases, e.g., the spectacu- 
lar merging system in the z=4.7 QSO BR1202-0725. Note that 
these SFRs are estimated under the assumption that the sub- 
mm sources are not lensed. This intense star-formation scenario 
is sup ported by the detection of the essential molecular r n aterial 



pp 

(e.g.. [Solomon & Vanden Boutl 120051: iMaiolino et al.L 12007 



Coppin et al. , 20081) . and the presence of luminous Polycy clic 



in these prodigiously star-forming galaxies. Assuming accre- 
tion onto a black hole at the Eddington rate, with efficiency 
e ~ 0.1, the typical observed z~2 optically-bright (Z?-band ab- 
solute magnitude — 27.5) QSO is hosted by a spheroid of mass 
M,ph ~ 10"MO. 

Previous redshifts of QSOs in the z~2 sub-mm de- 
tected sample were based on broad (5,000-10,000 km s"') 
high-ionization rest-frame UV lines that are observed in the 
optical, e.g., CIV, CIII], Lya, and NIV. Diff'erent broad 
emission lines in an ensemble of QSOs show systemati- 
cally different redshifts, with the high-ionization lines (e.g., 
CIV ^1549, CIII]/11909) displaying mean velocity diff^er- 
ences of ~500 km s~ ' and disp ersions (I cr) about this mean 
of 100-200 km s ' jTvtler & Fan! il992t IVanden Berk et al ' 



200 Ih . In addition it is not clear if the broad lines trace the 
systemic redshift of the host galaxy (as repr esented by the stel- 
lar content of the galaxy). jPenston ( 1977 ) proposed that the 
gas which emits narrow forbidden lines (e.g., [OIII]) in ac- 
tive galaxies is likely to give a good estimate of the systemic 
center-of-mass redshift. This low density gas comes from a re- 
gion ~ 1 kpc in extent and its motion should be dominated by 
the stellar density field in contrast to the broad line region gas 
which lies within a few parsecs of the central black hole. In lo- 
cal active galaxies, redshifts from these narrow (<1000 km s"') 
forbidden lines have shown agreement to <100 km s ' of the 
accepted systemic rest-frame determined from stellar absorp- 
tion features and HI 21 cm emission in host galaxies. 

Broad permit t ed lin es in QS Os, both rest-fra me UV lines 
(iBrotiierton et all 1 1994 and U/3 dHu et all l2008). can often be 
decomposed into a very broad component (>5000 km s ') and 
an intermediate-width component (~1500-3000 km s"'; typ- 
ically 0.4 times the width of the very broad component). The 
corresponding intermediate line region (ILR) is interpreted to 
be the outer regions of the broad line region (BLR) and may 
be dominated by infall dHu et all 120081) . The velocity shift be- 
tween the ILR emission lines and the [O III] line s is typically 
within a few hundreds of km s~^ (Huetal, 2008j). The former 
is thus a better tracer of the QSO systemic velocity in compar- 
ison to rest-frame UV lines. 

Black hole masses in QSOs can be estimated via 
scaling relations derived from reverberation 



scaling relations d e rived irom reverberation mapping 
dBlandford & McKeel 1 19821; IVestergaard & Petersonl I2OO6I) 
results: specifically via a single epoch measu rement of a 



aromatic hydrocarbon (PAH) emission ( Lutz et all 120081) in continuum luminosity (representing BLR size; iKaspi et al 



several of the sub-mm detected QSOs. The correlation between l200m and the FWHM of a broad permitted li ne, e.g., Hyg . 



PAH emission and sub-mm emission dLutz et all 120071 120081) 
supports star formation related heating of the dust rather than 
AGN h eating. The studies of optically bright, radi o-quiet QSOs 
at z~2 dPriddev et alll2003aHOmont et alll2003h were the first 
systematic studies of the star-formation propert ies of QSOs at 
the epoch that spans th e peak of AGN activity dSchmidt et al 



Ha, C IV (representing gas velocity in the BLR; iBentz et al 
20091) . The most direct and best characterized single-epoch 



method, derived from reverberation mapping results in nearby 
(relatively low luminosity) Seyfert galaxies, is via the pair 
Lsioo and FWHMh/j. The single-epoch black hole mass 
estimation base d on this pair is ex pected to be accurate to a 



tne epocn tnat sij ans tne peak or auin activity (lacnmidt et ai.L estimation base d on tnis pai r is expected to be accurate to a 
199U IWaU etalL l2005 b and resulted in the mm and sub-mm factor ~2, (e.g., iBentz et alj . l2009l) . The extrapolation of this 



detection of 18 QSOs in the redshift range 1.5 < z < 3.0. The 
observing approach adopted in the SCUBA (850jum) survey 
was to identify the very brightest (850//m flux >7mJy) of the 
sub-mm-bright QSO host galaxies with the aim of establishing 
a small sample suitable for follow-up observations in the pre- 
ALMA era to search for the molecular gas reservoir expected 



black hole mass estimator to high luminosity QSOs has been 
questioned on two main fronts: (a) is the relationship between 
continuum lumino sity and BLR size (the 7?blr -i- relationship; 
Kaspi et all 120051) unchan ged at higher lum inosities? Initial 



results with the C IV line (iKaspi et all 12007,) imply that it is 



relatively unchanged though a firm result using the HyS line 
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Standard Stars Observed 



Name 


R.A. 


Dec. 


Vhdio 


Spectral 


J,H,K 


Dale 




(J2000) 


(J2000) 


(kms"') 


Type 


(mag) 


Obsei'ved 


HD90508 


10:28:03.9 


+48:47:05 


-6.6 


F9V 


5.195,4.895,4.874 


March 2002 


HD93521 


10:48:23.5 


+37:34:13 


-16 


09Vp 


7.499, 7.647, 7.696 


March 2002 


HD 186980 


19:46:15.9 


+32:06:58 


+4 


07.5 III 


7.237, 7.272, 7.291 


March 2002 


1E0036. 1+3309 


00:38:47.8 


+33:25:43 




G2V 


9.057, 8.739, 8.607 


September 2002 



Note. — Columns are: f 1) star name; (2) and (3) position in J2000 coordinates (from SIMBAD); (4) heUocentric 
velocity (from SIMBAD); (5) spectral type; (6) I, H, and K magnitudes from the 2MASS point source catalog; (7) 
observation period at the TNG. 



is awaited ( Botti et al. . 2010*). The results in this work are 
relatively unaffected by small changes in the relationship - 
cur rently tested to con tinuum luminosities up to Log (A L,i) = 



and the Hamburg Bright Quasar Survey (HBQS; 'Hagen et al.l 



46 ( Bentz et al. , 20091) - as the QSOs discussed here have con- 
tinuum luminosities of Log (A L^) ^ 46.3^6.8; (b) should the 
FWHM of the H/3 line be estimated only from the Very Broad 
Line Region (VBLR) profile or from the combined VBLR and 
ILR profile? It is not uncommon that reverberation-mapped 
Seyferts show complex broad line profiles which cannot be fit 
with a single broad Gaussian. Further, complex gas motions 
including inflows an d outflows are seen in these Seyferts (e.g., 
Dennev et al. l2009l) . Even in these complex cases it appears 
that using a simplistic FWHM for the full H/3 profile (after 
subtraction of only the NLR emission) does not appreciably 
chang e the scahng r elatio ns used to estimate the black hole 



mass (iDennev et al 



,|2010^. On the other hand, there is some 
evidence that only a few local Seyferts (2 of 10 studied) show 
a true ILR component; for these, the inclusion of the ILR 
component in the FWHM measurement leads to erroneous 
black hole estimations (Zhu et al., 2009; Zhu & Zhang, .2010). 

An initial motivation of this work was the derivation of an 
accurate QSO redshift in order to obtain molecular gas detec- 
tions on instruments with bandwidth 1-2 GHz. This required 
redshifts accurate to, e.g., 1% at the 115 GHz CO J:l-0 Une. 
The recent availability of extremely wide-band receivers, e.g., 
Z-SPEC, ZEUS, GBT-Spec, has now relaxed this requirement 
for a very accurate redshift. Rest-frame optical spectroscopy 
of these unique sub-mm-bright QSOs also permits constraints 
on their black hole masses and accretion rates, thus allowing 
comparisons between galaxy and black hole growth. Further, 
an accurate systemic velocity allows better constraints on the 
kinematics of the nuclear gas. 

In this work we present near-IR spectroscopy probing the 
Ha line and nearby continuum in ten of the 18 QSOs in a z~2 
sub-mm detected QSO sample. The sample selection, observa- 
tions, and data processing are summarized in Sect. |2] results 
are presented in Sect. |3] and a brief discussion and conclud- 
ing remarks are presented in Sect. Luminosity distances of 
the QSOs were calculated using a Hubble constant Ho= 72 
km S-' Mpc-', = 0, Qn, = 0.3, Qa = 0.7, and qO=-0.5. 

2. Sample, Observations, and Data Processing 



The sample is d r awn from iPriddev et al.l (l2003al) and 



Omont et all (l2003h . IPriddev et al.l (l2003a') surveved 57 opti 



caUy luminous (Mb < -27.5), radio-quie t QSOs at z~2 from 
flie Large Bright Quasar Survey (LBQS; iHewett et al.i Il995l) 



19991: lEngels et al.L Il998t iReimers et al.L 11995) at 850;um with 
JCMT/SCUBA. Here radio quiet was defined by a non- 
detection i n the NRAO VLA Sky Survey (NVSS), i.e. Si.4ghz < 
1.5 mJjOl IPriddev et al] (l2003al) detected 9 QSOs at a sig- 
nificance of 3cr or greater; these have SSQfim fluxes in the 
range 7-17 mJy. The remainin g QSOs have 3cr fl ux density 
limits in the range 6-10 mJy. lOmont et al.l (l2003l) surveyed 
35 optically luminous (Mb < -27.0) radio-quiet QSOs with 
1.8 < z < 2.8 at 1.2 mm with IRAM/MAMBO and detected 
9 QSOs with 1.2 mm fluxes of 3.2-10.7 mJy; the remaining 
QSOs have 3 i t flux density limits in the range 1.8^ mJy. The 
Omont et al. (2003) sample was drawn from the AGN cata- 
log of Vero n-Cettv & Veron (.2000) and included a variety of 
QSO surveys, among t hem the HBQS, t he LB QS, the Second 
Byurakan Survey (SBS; Stepanian et al., '2001') and the Hewitt- 
Burbidge QSO compilation (HB89; Hewitt & Burbid^ e. 198 9). 
Of the 83 unique QSOs surveyed by Priddey et al. and Omont 
et al., nine QSOs were observed by both. Of these nine, eight 
were not detected in both studies while HS B 1049+4033 was 
detected by Omont et al. but not detected by Priddey et al. 
Thus, combining the samples of Priddey et al. and Omont et 
al. resulted in 83 observed QSOs of which 18 were detected 
at 850/im or 1.2 mm. Given the relatively high upper-limits 
of the QSOs not detected at 850yum and/or 1.2 mm (in many 
cases the 3cr flux upper-limits are higher than the fluxes of 
the weaker sub-mm detected QSOs), we refrain from using 
the term sub-mm bright, and instead refer to the 18 QSOs as 
sub-mm detected QSOs, and the remaining as sub-mm non- 
detected QSOs. We have obtained near-IR spectra of ten of 
these 18 sub-mm detected QSOs, listed in Table|2] which there- 
fore form the sub-mm detected QSO sample discussed in this 
work. 

Observations were conducted on 2002 March 28, and 
2002 September 20 at the 3.6 meter Telescopio Nazionale 
Galileo (TNG) at La Palma, using the Near-IR Camera and 
Spectrograph (NICS) instrument with a detector b ased on a 
HgCdTe Hawaii 1024x1024 amy (iBafl^a et al.Ll200ll) . On 2002 
March 28, we used TNG/NICS wifli the Kb grism and a slit- 
width of 0'.'7 or 1". This yielded the wavelength range 1.91jum 
to 2.38yum at a dispersion of 4. 5 A per pixel i.e. R~1250 for a 1" 
slit. On 2002 September 20, we used TNG/NICS with the HK 
grism and a slit of 1 " width. This yielded the wavelength range 
1.37yum to 2.49//m at a dispersion of llA per pixel. Standard 

' The assumption of a spectral index a = -0.5 or steeper would 
imply SgsoGHz < 0. 1 mJy. 
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calibration frames were taken at the beginning and end of the 
observing night. The QSOs were observed using typical inte- 
gration times of 120 sec in an "ABBA" nodding pattern. The 
distance between "A" and "B" positions was roughly 15", and 
was varied slightly on each nod. Total integrations were typi- 
cally 30^0 min on-source per QSO. At least two telluric 'stan- 
dard star' observations were obtained on each night (Table [T}. 

All spectra were first "straightened' using 54 strong night 
sky lines. Xe and Ar lamp exposures, taken at the beginning and 
end of the night, were used for initial wavelength calibration, 
with individual lines identified using the line-lists (in vacuum 
wavelengths) available from the ESO/ISAAC websit^ The Ar 
spectrum is richer in lines in the case of the Kb grism: of these 
we used 1 1 reliable lines. The Xe spectrum contributed three 
additional lines to the Kb grism observations. These 14 lines 
ranged between 19823 A and 23200A, which nicely covered 
the full wavelength range of our setup: ~19157A to 23800A. 
Two QSOs - [HB89]0933-h733 and SBS B1408-H567 - have 
Ha lines which fell at the red-most end of the CCD: in these 
cases we checked that using additional weaker Ar and Xe lines 
in the red-most region gave a similar wavelength calibration 
solution. The final wavelength solution from the lamps had an 
rm.s. about 0.3A, with almost all points well within 0.5A of 
the fit. In the case of the HK grism, we used only the Ar lamp, 
which yielded 34 usable lines, between 13179A and 23973 A. 
The final fit had an r.m.s. of 4A; at the blue end the points were 
all within ±5A while at the red end the points were all within 
+ 10A. The wavelength calibration was checked both with the 
calibration star spectra and with night sky lines: in all cases the 
observed lines were within 1 A (March 2002) or 3 A (September 
2002) of their expected values. In the case of data taken in 
March 2002, the expected wavelength of the night sky lines 
was shifted by a consistent amount (within 3A) in observations 
made towards the end of the night. The latter shift was probably 
due to a slight change in tilt of the grism with telescope move- 
ment, and we corrected the wavelength solution of the target 
objects with an equivalent zero-point shift. The accuracy of the 
wavelength calibration was therefore within 2A or ~30 km s"' 
at the typical observed wavelengths of the Ha line. 

Finally, QSO spectra were continuum subtracted assuming 
a linear baseline and the Specfit dKrissl 1 19941) package within 
IRAF was used to obtain Gaussian and Lorentzian fits to the 
Ha line profile. The very broad profiles of the lines results in 
an uncertainty of about 5-1 OA in the fitted central wavelengths 
(10-15A for the September 2002 data due to its lower spectral 
resolution). This is the dominant source of measurement error 
in the redshift determinations. 

We note that NICS data was taken soon after instrument 
commissioning, at a time when telescope and instrument point- 
ing was still being fine tuned. The telescope nodding along 
the slit was therefore slightly inaccurate, and the counts on the 
standard stars and QSOs varied, in the worst case, by up to 20% 
between exposures. Flux calibration was performed by compar- 
ing the best spectra of the standard stars to the expected spectra 
as estimated from black body fits to the V- to K-band mag- 
nitudes of the standard stars. The flux calibrated QSO spectra 



were then used to estimate the J, H, and K-band magnitudes of 
the QSOs. We found that the spectral fluxes had to be increased 
by ~30% in order to be consistent with the observed 2MASS 
magnitudes of the QSOs. This correction, likely the result of 

slit losses, was thus applied to all QSO spectra. 

We use the primary Q SO samples of 'Shemmer et al.l 



(I2004D . iNetzer et all (120071) and Alexander et al. (2008) to 
compare and contras t the properties of our QSO sample. 
Shemmer et al. 1 (12004 ') present the results of new rest-frame op- 
tical spectra of a sample of 29 high-reds hift ( ty pically z~2.3 ) 



http://www.eso.org/instruments/isaac/tools 



luminous (L > 10"*^ erg s ') QSOs and iNetzer et all (120071) 
present the results of new rest-frame optical spectroscopy of 15 
luminous QSOs at z~2. 3-3.4. Alexander et al. (2008) present a 
sample of z~2 sub-millimeter detected galaxies (SMGs) which 
exhibit broad Ha or broad HjS emission. All three works use 
the results of the rest-frame optical spectroscopy to derive es- 
timates of the black hole mass and related quantities which are 
directly comparable to our results. In all cases we convert listed 
fluxes into luminosities using the cosmologi cal p a ramete rs of 
this work. We also use the results of ICoppin et al.l (I2OO8I) who 
present rest-frame optical spectroscopy of a sample of sub-mm 
detected QSOs, four of which are also in our sample. 

3. Results 



3.1. Ha Profile Fits 

Spectra covering the region near the Ha line are shown in Fig.[T] 
for all observed QSOs. None of the sample QSOs have Ha pro- 
files which can be satisfactorily fit with a single Gaussian. Two 
components, a very broad (4500-9000 km s"') Gaussian and 
an intermediate-width (typically 1500-3200 km s"') Gaussian, 
are instead required. Narrow lines (<1000 km s"') are not 
clearly detected in any of the sample QSOs, though there is 
a very tentative detection of narrow Ha + [N II] lines in LBQS 
B0018-0220 (Fig.Q. Both the very broad and the intermediate- 
width Gaussian components are centered on more or less the 
same wavelength (diff'erences of <10 A or <120km s"'), for all 
QSOs except SBS B1542+541, which has exceptionally wide 
ILR and VBLR Ha components and an obvious blue shoulder, 
and HS B0035-I-4405, for which a reliable double-Gaussian fit 
was difficult given that the red wing of the Ha profile extended 
beyond the observed range. In the two latter cases, the differ- 
ence between the fitted central wavelengths of the very broad 
and intermediate-width Gaussian components, 40A, represents 
a velocity difference of about 500 km s ' . On the other hand, a 
single Lorentzian profile provides a satisfactory fit to the bulk 
of the Ha line profile in most QSOs, with the main drawback 
of potentially overestimating the line profile in the outer wings. 
This single Lorentzian fit to the Ha profile has a central wave- 
length similar to, or the same as, that of the intermediate-width 
Gaussian component in the double-Gaussian fit. 

We calculate the systemic redshift of the QSO from the cen- 
tral wavelength of the intermediate-width (ILR) component of 
the double-Gaussian fit to Ha. The total flux of the Ha line 
is taken as the summed flux of the double-Gaussian fit. Given 
our non-detection of any narrow (NLR) component, the mea- 
sured Ha fluxes and widths are dominated by the BLR (VBLR 



Gustavo Orellana et al.: Sub-millimeter Detected z ~2 Radio-quiet QSOs 



5 




-10 -5 5 10 15 -10 -5 5 10 15 -10 -5 5 10 15 -10 -5 5 10 15 -10 -5 5 10 15 



Velocity x 1000 km s"' 

Fig. 1. Near-IR spectra of the sample QSOs centered on the Ha emission line, along with our double-Gaussian fits to the line 
profile. The two Gaussians represent an intermediate-width component (corresponding to the ILR), and a very broad component 
(corresponding to the VBLR; see text). For each QSO the upper panel shows the normalized QSO spectrum (black) along with 
the double-Gaussian fit (grey); the lower panel shows the residual spectrum after subtraction of the double-Gaussian fit. The zero 
velocity is based on the central wavelength of the intermediate- width component of the double-Gaussian fit (see Table |2]i. The 
ticks below each QSO name indicate the expected locations of the [N II] doublet. 



plus ILR). The estimation of black hole masses for the sam- 
ple QSOs requires the measurement of the full-width at half 
maximum (FWHM) of the broad component of the Ha emis- 
sion profile, i.e. the FWHM of the Ha profile after subtrac- 
tion of any narrow (<1000 km s ') component. Given that a 
single Gaussian does not provide a good fit to the line profile 
and that the narrow component is not detected in any QSO, we 
calculate the FWHM of the Ha line profile from the summed 
double-Gaussian fit, and refer to this quantity as FWHMea. 
This FWHMh„ is within 50-100 km s ' of the FWHM of the 
single Lorentzian fit for all QSOs except SBS B1542+541 for 
which the difference is 620 km s ' (Table|2]l. 

The double-Gaussian fits to the Ha lines are shown in 
Fig.[T]for the sample QSOs, and the corresponding line centers, 
FWHMs and fluxes of each Gaussian component are listed in 
Table 12] Values of the FWHM of the single-Lorentzian fit and 
of FWHMh^ are also listed in the same table. 



3.2. Systemic Redshifts and Velocity Offsets 

In this section we first compare our redshift estimations to pre- 
viously obtained redshifts based on rest-frame-optical spec- 
troscopy. We then analyze differences between redshifts de- 



rived from Ha and from rest-frame UV spectroscopy for the 
sample QSOs. We further study in detail the velocity offsets 
between individual emission lines for the subset of the sam- 
ple QSOs which have high quality rest-frame UV spectra. 
Finally, we contrast the sub-mm detected and sub-mm non- 
det ected QSOs in the p arent samples of Priddev et al. (2002i) 
and lOmont et al.l ( 2003h with respect to velocity offsets seen in 
different emission lines. 

Four of our sa mple QSOs are also in the sample of 
Coppin et al.l OOOSh . who present results of near-IR spectro- 
scopic observations which include the [OIII] /i5007, H/3 and 
Ha emission fines, and observatio ns of the J:3-2 or J:2-l tran- 
sition of CO. ICoppin et al.l ( l2008h derived QSO systemic red- 
shifts via fits to the HjS and [O III] fines. A comparison between 
our Ha determined redshifts and those of ICoppin et al.l (l2008b 
shows an agreement of within ~300 km s~' in two QSOs, but 
larger diff'erences in HS B1611+4719 and SBS B1542H-541 
(Fig. IS. The latter is the QSO in which both ILR and VBLR 
Ha components have exceptionally large widths and in which 
a notable Ha blue shoulder is seen. Our Ha derived red- 
shift is therefore probabl y not as accurate as the [OIII] and 
CO derived redshifts of (ICoppin et al.l 120081) . In the case of 
HS B1611H-4719, our relatively refiable fit to the Ha line re- 
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suits in a Ha redshift which is offset by 900 km s ' from the 
Hj6 and [OIII] derived redshi ft, but offset by only 280 km s ' 
from the CO derived redshift feoppi n et aObOOSl) . 

Our Ha derived redshifts are systematically higher (Fig. |2) 
than those derived from optical spectroscopy (i.e. based on 
rest-frame UV lines). The rest-frame UV redshifts originate 
from very heterogenous observations with diverse resolutions 
and spectral ranges: these are taken from, in order of prefer - 
ence: SDSS, Hewett et dj (Il995h . iHewitt & Burbidgd Jl989l) . 
IStepanian et al. ( 200 1 ), and the Ha r nburg Bright Quasar Survey 



( iHagen et all Il999t lEngels et all Il998l) . Of these, only the 



SDSS redshifts have quoted uncertainties; the other sources do 
not explicitly specify uncertainties. In the latter cases we use 
indicative uncertainties estimated from the spectral resolution 
and/or the number of significant digits quoted. The median dif- 
ference is Az = 0.018 and the mean difference is < Az>- 0.021, 
corresponding to AV ~ 6300 km s ' at wavelengths near Ha. 
The largest redshift differences are typically for Ziestuv from 
the Hamburg Bright Quasar Survey which has a low resolution 
(typically 15A), limited (and relatively blue) spectral coverage, 
and few significant figures in the reported redshift (with no er- 
rors specified), while the smallest differences are typically for 
Zrestuv ffom the SDSS. Thus the redshift differences larger than 
Az ~ 0.02 are likely the result of the low spectral resolution, the 
limited set of lines used, and the low precision, of the redshifts 
derived from rest-frame UV lines. Limiting ourselves to the 
QSOs with accurate rest-frame UV redshifts the typical rest- 
frame optical to rest-frame UV redshift difference is ~0.013, 
or~4000kms-i. 
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Fig. 2. Comparison of redshifts derived from the Ha line and 
those derived from rest-frame UV spectroscopy (SDSS as tri- 
angles and HBQS as diamonds). The median difference is Az = 
0.018. Error bars in y are indicative for non-SDSS QSOs (see 
text); error bars in x are roughly the size of the symbols. 



In order to better determine the velocity offsets of individ- 
ual rest-frame UV lines we limited ourselves to the sample 
QSOs for which high quality rest-frame UV spectra exist. To 



this end, we searched the Sloan Digital Sky Survey (SDSS; 
Data Relea se 7; DR7) archive fo r sp ectra of all QSOs in the 
samples of'Priddev et al.' ('2003a') and bmont et al.1 (l2003h and 
also those of Coppin et al. (2008). We use the DR7 spectra and 
their associated "SpecLine" tables which tabulate the SDSS au- 
tomated fits to all detected emission and absorption lines in 
the spectra. We selected emission lines from the SpecLine ta- 
bles using the following criteria: (a) line width (cr) between 
20A and 200A. The lower limit avoids spurious line detections 
and the upper limit avoids lines which are suspect and/or can- 
not yield sufficiently accurate line centers; (b) error in /lobs less 
than the line width (cr); (c) error in /lobs less than 0.5% of /lobs- 
Since the SDSS line identification and fitting are automatic, 
there are some errors in the database of line parameters. We 
therefore continuum subtracted the SDSS spectra using a poly- 
nomial baseline fit and revised all lines selected above, man- 
ually re-fitting the line centers and widths if necessary. This 
led to a change in the central wavelengths and widths of some 
lines and the elimination of other lines. The re-fitted values of 
/lob.s and their Icr eiTors were used to derive relative velocities 
between emission lines. When a Ha or [OIII] derived redshift 
was available, this was used as the systemic redshift; else only 
relative velocities between rest-frame UV lines were analyzed. 
A comparison of the velocity offsets betw een emission 



lines for the QSOs in our sample and those of ICoppin et al 
(120081) are shown in Fig. |3] Almost all rest-frame UV 



emission lines show large blue-shifted velocities; several 
factors higher than typically seen in QSOs (Tvtle r & Fan! 



1992t IVanden Berk et al.L 120011) . with the most extreme blue 



shifted velocities seen in HS B 1049-1-4033. We note that 
HS B 1049+4033 is not extreme in i ts sub-mm p r opertie s: it 
was detected at 1.2 mm (3.2 mJy; lOmont et al.L 120 03) but 
not detected at 85 0yum (3.9+3.2 mJv; IPriddev et all i2003al) . 



Vanden Berk et al.l (POOl) present median velocities for emis- 



sion lines relative to [OIII] A5001 for over 2000 QSOs from 
the SDSS over the redshift range 0.04 < z < 4.8. Relative 
to the [OIII] line, they find median velocity differences of 
27 km s-i for Ha, -70 km s ' for NV, -563 km s~i for 
CIV, -224 km s ' for CIII], and 161 km s ' for Mgll. 
Relative to Ha (or in the absence of Ha, relative to [OIII]) 
derived redshifts, the sub-mm detected QSOs in Fig. |3]have 
CIV velocities of -1230 km s"' to -3420 km s"', and CIII] 
velocities between -900 km s"' to -2070 km s"' except 
for SMMJ 123716.01+620323.3 which has a CIII] velocity 
of 1500 km s~' relativ e to the [QUI] systemic velocity pub- 



lished by ICoppin et alj (120081) . Note that the sample QSOs 
show small blue-shifted velocities for Mg II and Lya instead 
of the small median redshifted v elocities seen in the QSO sam- 
ple of IVanden Berk et"aL (2001); the Lya velocity offsets seen 
in Fig. [3] are not reliable as the Lya line center is difficult to 
ascertain given its blue-ward absorption and the contamination 
from the N V line. 

In the case of the parent samples of sub-mm detected 



QSOs vs. sub-mm non-detected QSOs (iPriddev et al.L l2003a : 



Omont et all l2003h . 40 of the 89 Q SOs have SDSS spectra 



and also appear in the compilation of lShen et al.l ((2010). In this 
larger sample, we do not have reliable QSO systemic veloci- 
ties via the Ha or [O III] line for all QSOs and therefore only 
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Fig. 3. Velocity offsets of rest-frame UV, rest-frame optical, and CO emis sion lines for sub-mm detected QSOs. The QSOs plotted 
are present in at least two of the following: our sample, the ICoppin et al. (2008) sample, and the SDSS spectroscopic database. 
The lines are ordered left to right in increasing rest wavelength, individual lines are labeled, and different symbols, as indicated 
on the figure, are used for each QSO. Only lines with reliable fits to their line centers are included (see text). The QSO systemic 
velo city is assumed to be that of the intermediate- width (ILR) component of the Ha line if available, or else that of the [O III] line 
from Coppin et alj ( 2008h. Error bars have been displaced slightly along the x axis for better visibility. The QSO SBS B 1542-1-541 
is in our sample. 'Copp in et all (l2008h . and has an SDSS spectrum. Two QSOs, HS B1049-H4033 and HS Bl 103+6416, are both 
in our sample and have SDSS spectra. Two QSOs, HS B1002-H4400 and SMM J123716.01+620323.3, are both in Coppin etal] 
(l2008h and have SDSS spectra. Three addit ional QSOs, SBS B1408-H567, HS B161 1+4719, and J164914.9+530316, are both in 
our sample and that of ICoppin et al.l (l2008b . 
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Fig. 4. A comparison of the distributions of Left: CIV velocity relative to the SDSS-derived systemic velocity of the QSO, 
Middle: CIV velocity relative to the Mgll velocity, and Right: black hole mass, contrast i ng the sub- mm detected QSOs (solid 
lines) and sub-mm non-detected QSOs (dashed lines) from the samples of iPriddev et al.l (12003 a) and O mont et al. ( 2003h . The 
histograms have been shifted slightly in the x-axis for legibility. Values used are from the comp ilation of | Shen et al. [(120101) except 
for the inclusion of the 10 Hcc-based black hole masses presented in this work. In contrast to lShen et al. I (l2010h we use positive 
velocities to indicate redshift. 



tested for relative velocity offsets between rest-frame UV lines. 
No significant difference between the sub-mm detected and 
sub-mm non-detected sub-samples was found. We also used 
the ve locity offse t s of th e C IV and Mg II emission lines listed 
in the Shen et al. ( 2010h compilation of QSO properties (these 
are the only two rest-frame UV emission lines for which ve- 
locity offsets from the SDSS-derived systemic QSO ve locity 
are listed in the compilation). Note that lShen et al. I (l2010l) used 



negative velocities to indicate redshift, while we use positive 
velocities to indicate redshift. The velocity offsets of the C IV 
lines are not significantly different for the sub-mm detected 
and non-detected subsamples (left panel of Fig. HJi. Velocity 
offsets for the C IV and Mg II lines were calculated relative to 
the SDSS-derived redshift of the QSO. This SDSS-derived red- 
shift - based on a cross correlation of narrow emission and ab- 
sorption lines - does not necessarily trace the systemic redshift 
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of the QSO as it is likely biased by the large velocity offsets 
in rest-frame UV emission-lines as discussed previously. The 
Mgll line typically shows velocity offset s which are sni aller 
than those of the C IV line, e.g., Vanden B erk et al.l ( 1200 lb and 
Fig. [3] We thus also compare the velocities of the C IV lines rel- 
ative to those of the Mg II lines, which potentially better repre- 
sent the velocity offset of the C IV lines from the true systemic 
QSO redshifts (middle panel of Fig.|4|i. While both subsamples 
show large blue-shifted velocities, there is no clear difference 
between the two subsamples. 

3.3. Black Hole Masses and Related Properties 

The relationships between the three primary measurements 
used in this section - Ha luminosity, FWHMHa, and continuum 
luminosity - are shown in Fig.|5] The Ha and continuum lumi- 
nosities are proportional (left panel of Fig.|5), with a continuum 
luminosity a factor of ~2 higher than would be expected from 
the extrapolation of the same relationship app lied to a large 
samp le of SDSS galaxies (log L^m ~43-45; lOreene & Hoi 
20051). This correlation is unlikely to be due to a distance bias: 
the QSOs in the sample have similar redshifts, and a simi- 
lar correlation is seen in the equivalent flux-flux plot. There 
is no clear relationship between FWHMh» and both the lumi- 
nosity of the Ha line and the continuum luminosity (Fig. |5]l. 
Compared to the s ub-mm detected z~2 gal axies with broad Ha 
or broad H/3 from [Alexander et al.l (l2008h . our sample QSOs 
all have higher Ha luminosities, but show a considerable over- 
lap in Ha line widths. In continuum luminosity and FWHAIho, 
our sample QSOs are very similar to the QSO sample stud- 
ied by Shemmer et al. (2004), but typically an order of magni- 
tude more luminou s in the continuum than the QSO sample of 
Netzeret al.1(l2007h . 

Estimates for the black hole masses, Eddington frac- 
tions (i.e. the fractional Eddington rate), mass accretion 
rates, and FIR luminosities are listed in Table |3] for the 
10 sample QSOs. We estimated the black hole masses and 
Eddington fractions of the sample QSOs from the empir- 
ical scalings derived from reverberation mapping results. 
Specifically, we used eqn. 1 of Greene & Ho (2005) to con- 
vert FWHMHff to FWHMh/3 and then estimated the black 
hole mass from FWHMh/j and L5100 following eqn. 5 of 
Vestergaajd & Peter son (2006 [). Following K aspi et al. (2000), 
Marconi et all (l2004l) . and iNetzer & TrakhtenbrotI (l2007l) we 
use L5100 - ^ L,i with A = 5100A, and LboI - 1 x Lsioo- 
Eddington fractions were computed assuming this bolometric 
luminosity, and the mass accretion rate was computed assum- 
ing a 10% accretion efficiency (e = 0.1). We calculated the FIR 
luminosity from the sub-mm flux (850/vm or 1 .2 mm) following 
McMahon etld] (fl999l) and using p - 1.5 and dust temperature 
Td = 40 K. Finally, the star formation rate ( SFR) was calculated 
from the FIR lumino sity following McMah on et al.l 61999.) and 
using the iKennicuttId 19981) relation, i.e. a - 1.72. 

The relations we used are therefore: 



FWHMh^ = (1.07 + 0.07) X 10^ 



— ; r km s 

103 km s- 



Mbh/Ms = 10' 



6.91 



FWHMug 



,1000kms-i 

Ifidd - LBol/LEdd = 7 Lsioo/LEdd 



L5 



\0.5 



100 



1044 erg s- 



dt e\ 1046 erg s^V \yr-^ 



= '''' ^ l353 GHi) 



3.3 X 



Ma 



\ 108 Mo 



X lO'^L, 



SFR = 1.72 X 10 



-10 



^FIR -\ 

—Moyr 



The lower rest-wavelength limit of our near-IR spectra is 
typically 5500A to 6000A. Given that the spectral slope in the 
5500A-6500A range is relatively flat for the QSOs in our sam- 
ple we simply approximate that the 5100A flux is the same as 
the 5500A or 6000A flux. FWHMh„, the FWHM of the Ha 
line profile as measured from the sum of the double-Gaussian 
fit, was used in the black hole mass related calculations. If the 
FWHM of only the very broad (VBLR) Gaussian component 
of the Ha line is used instead, the resulting black hole masses 
would increase by a factor of ~5. We also estim ated the blac k 



hole masses following the scaling relations in IMarconi et al 



(2 0081) with a nd wi thout radiation pressure effects, and in 
iGreene & Ho (2005); in all cases except one, the black hole 
mass estimates remain within a factor of 2-3 for almost all 



QSOs in the sample. The exception is the lMarconi et al.l (l2008l) 
relation which includes radiation pressure effects: here the high 
continuum luminosities of the sample QSOs result in the radi- 
ation term dominating the black hole mass estimate, i.e. the 
estimate is relatively independent of line width. The use of 
FWHMho a nd Lhq to estimate t he black hole mass (follow- 
ing eqn. 6 of lGreene & H o[ l2005l) leads to slightly lower black 
hole masses in all QSOs (Table O due to the slightly higher 
than usual luminosity of Lsioo compared to Lua (Fig.jSj. 

Errors in the black hole mass were estimated both by Monte 
Carlo simulations assuming Gaussian errors for the input quan- 
tities and by the analytical cal culation of the error propagation 
(iBevington & Robinsonlll992 ). Both methods resulted in con- 
sistent error estimates, with errors around 10%. Note, however, 
that these errors represent only the measurement errors and the 
true err or in the black hole mass is h kely closer to a factor 
4 (e.g., IVester^aard & Petersonl l2006l) . For the other derived 
quantities (e.g., Indd) errors were calculated by analytical error 
propagation. 

Two of the sample QSOs have bl ack hole es timates de- 
rived from the Mgll and CIV lines in lShen et al.1 (I2OIO). For 
HS B 1049-1-4033, our estimated value lies close to the values 
(9.5 to 9.8 with typical error of 0.07 as estimated from three dif- 
ferent scaling relationships) estimated from Mg II and signifi- 
cantly lower than the value (10.21+0.06) estimated from C IV. 
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Fig. 5. The relationship between continuum luminosity (A at 5100A; see text), and the width and luminosity of the broad 
component (ILR+VBLR) of the Ha line. QSOs in our sample are plotted with filled circles and their corresponding Icr error 
bars. The line in the left panel corresponds to the relationshi p between these quantit ies found by Greene & Ho (2005) in a sample 
of SDSS galaxies (see text). In the middle panel data from lAlexander et al. I (l2008h on the sub-mm detected (SMG) galaxies in 
the Chandra Deep Field North are shown as triangles (obscured ULIRGs with broad Paa) and squares (SMGs with broad Ha 
or broad H/3). The right panel includes the QSO samples of Shemmeret al. (2004, triangles) and Netzeret al. (2007, squares); 
quantities for these were recalculated using the scaling relations and cosmology used in this work, and smaller symbols are used 
for QSOs with redshifts >2.5. 



For HS B 1 103-H6 416, our est i mated value is significantly lower 
than that fisted in lShen et alJ (l2010l) : 10.4+0.06 and 10.3+0.03 
for Mg II and C IV, respectively. We remark that that using the 
FWHM of the very broad ( VBLR) Gaussian component to cal- 
culate the Ha based black hole mass leads to a better agreement 
with the C IV and Mg II based black hole mass estimates. 

The left panel of Fig. |6] contrasts the relationship between 
the fractional Eddington rate and black hole mass of sub-mm 
detected QSOs in our sample to two o ther QSO samples at z~2; 
those of IShemmer et all hoo4 and iNetzer et al] (|2007|) . AU 
three samples have a similar range of black hole masses, though 
the QSOs of the Shemmer et al. sample and our sample are ac- 
creting at relatively large Eddington fractions. Combining this 
result with that seen in the right panel of Fig. |5] shows that 
our sub-mm detected QS O sample is virtually in distinguishable 
from the QSO sample of IShemmer et al. I (l2004 . This result is 
not very sensitive to the exact scalings used in the black hole 
mass estimation (see the arrows in the figure): if the scaling 
relations are modified, all QSOs in the figure will shift their 
position in a relatively consistent manner Note that the scaling 
relations we use imply that L/LEdd °^ '^s foo (FWHM//(j,)"^ "^ and 
Mbh °^ L^'foo (FWHMffa)^ Thus Fig.|6]plots highly related 
quantities, so that any subtle correlations or anti-correlations 
seen are to be treated with caution. As mentioned in Sect. [T] 
it may be more correct to estimate QSO black hole masses 
from the FWHM of the very broad (VBLR) component in- 
stead of FWHMHff. The right panel of Fig. |6] demonstrates 
the implications of using the VBLR width in the black hole 
mass estimation: black hole masses increase, and fractional 
Eddington rates decrease, by a factor of ~5 . The compari- 
son QSOs dShemmer et all |2004 iNetzer et all l2007l) will also 



change positions in this plot, but in an unknown way as their 
VBLR widths are not listed by the corresponding authors. 

The relationship between FIR luminosity (proportional to 
the star-formation rate) and the mass accretion rate onto the 
black hol e (proportional to Lbqi) for the QSO sample is shown 
in Fig. |2] i Trakhtenbrot & Netzet (l2010l l) have studied a large 
sample of Seyferts and QSOs with redshifts 0.2 < z < 3, bolo- 
metric luminosities in the range 2 x 10"*^ - lO'*'^ erg s"', and 
SFRs in the range 10 ' - lO'* Mo yr ', and find that these 
can be well characterized by a single relationship (black line 
in Fig.|71l: 



SFRiM^yr-^) ^ 32.8 



LboI (erg/s) 
1046 



0.7 



Our sample QSOs are found close to the highest end of this 
relationship and off'set from it in having SFRs almost one mag- 
nitude hig her. The sample o f 12 millimeter-bright z~2 QSOs 
studied bv iLutz et al.l (120081) with Spitzer are also shown for 
compariso n : four of their targets are also present in our sample. 



Lutz et al.l (l2008h estimated SFR and bolometric luminosity 



from PAH-emission and 60;um continuum, respectively, lead- 
ing to some diff'erences from the values derived by us. 

Finally, we compare the black hole mass distributions of 
the parent sample of sub-mm detect ed and non-detected Q SOs. 
Of the 89 QSOs in the samples of IPriddev etaP (l2003al) and 



(I2010h . IShen et al 



lOmont et al.l ( 2003 ), 40 app ear in the compilation of lShen et al 



m 



(120101) list, among other quantities, the 



black hole masses derived from the C IV and Mg II lines and 
nearby continuum. The black hole mass distributions of the 
sub-mm detected and sub-mm non-detected QSOs (right panel 
of Fig. are not significantly different. To realize this figure, 
we derived the black hole mass of the QSO by using, in order 
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Fig. 6. Left: The relationship between Eddington fraction and black hole mass contrasting QSOs in our sub-m m sample (filled 
circles plotted with their Icr errors) to the QSO samples of .Shemmeret al. (2004. triangles) and iNetzer et alj (12007. squares). 
Quantities for the latter two samples were recalculated using the scaling relations and cosmology used in this work, and smaller 
symbols are used for QSOs with redshifts >2.5. The thin and thick arrows show the typical shift of the IShemmer et al.l (l2004h 
and .Netzer et al.l (12007 ) points, respectively, if the QSO black hole masses are calculated via the equations used by these authors. 
Right: As in the left panel, but only for QSOs in our sub-mm sample and in the case that only the very broad (VBLR) Hor 
component is used in the black hole mass estimation. 



of preference, our Ha derived black h ole mass, or the M g 11 or 
C IV derived black hole mass listed in Shen et al. I (l2010h . This 
resulted in black hole masses for 14 sub-mm detections and 33 
sub-mm non-detections. 



4. Discussion and Concluding Remarl<s 

We have observed the Ha line and neighboring continuum in 
lOofthe 18 sub-mm detected QSOs from a parent sample of 83 
optically-bright QSOs at z~2. The properties of these 10 QSOs 
have been compared with a sample of sub-mm detected galax- 
ies and two other samples of optically-bright z ~ 2 QSOs. We 
have also attempted to compare the 18 sub-mm detected QSOs 
with the 65 sub-mm non-detected QSOs in the parent sample. 
This comparison is limited by two factors; (a) the lack of ac- 
curate systemic velocities for the sub-mm non-detected QSOs 
and. (b) the relatively low sub-mm fluxes of some of the sub- 
mm detected QSOs coupled with the relatively high detection 
limits of the sub-mm surveys makes it difficult to clearly sep- 
arate sub-mm bright QSOs from sub-mm dim QSOs. That is, 
it would not be surprising if some of the QSOs not detected 
in the sub-mm have sub-mm fluxes similar to the weakest sub- 
mm detected QSOs. Nevertheless, the parent sample used here 
(18 detections in a sample of 83 QSOs) is sufficiently large to 
obtain meaningful comparisons between sub-mm-bright QSOs 
and other QSOs. 

The Hff line profiles can be decomposed into contribu- 
tions from boffi a very broad line region (VBLR) and an 
intermediate-width line region (ILR), similar to the cases of 
Hy6 and CIV in QSOs. Both components show similar red- 
shifts in most of the QSOs, and the ratios of the VBLR to ILR 
widths (1.7-5 with a median of 3.0) are not different from the 



typica l ratios (~2.5) found for HyS profiles in QSOs (iHu et al 
2008h . The ILR has been interpreted to be the outer portion of 



the BLR and may originate from a physically distinct region 
(IHu et al.l |2008) . Its inclusion in the black hole mass estima- 



tion has been questioned (IZhu et al.Ll2009h . Following previous 



works, we estimate the QSO black hole masses using the com- 
bined VBLR+ILR width, but also discuss the implications of 
using only the VBLR width. 

We have made the reasonable assumption that the 
intermediate-width (ILR) component of the broad Ha emission 
closely traces the true systemic velocity of the QSO. A blind 
comparison of Ha-based redshifts derived here with previous 
redshifts based on high-ionization rest-frame UV-lines shows 
larger than expected redshift differences (Fig. |2]i. However, 
the largest differences are seen in QSOs whose rest-frame UV 
redshifts are imprecise, inaccurate, and/or for which redshift 
errors have not been provided. For example, the HBQS and 
LBQS quote redshifts for these QSOs to the second decimal 
place (i.e. an accuracy of ~1500 km s"' if the redshift val- 
ues were rounded or ~3000 km s ' if truncated). Limiting 
ourselves to the QSOs with accurate rest-frame UV redshifts, 
the typical rest-frame optical to rest-frame UV redshift differ- 
ence is ~0.013, or ~4000 km s Further limiting ourselves to 
only sample QSOs with reliable rest-frame UV spectroscopy 
from SDSS, we find larger than average velocity offsets (up to 
3400 km s"' ) for individual rest-frame UV lines (Fig. |3)- Given 
the large obscuration expected in sub-mm QSOs (nuclear dust 
would obscure gas on the far side of the AGN), these veloc- 
ity offsets most likely trace ionized gas outflows. The lack of 
accurate redshifts for the sub-mm non-detected QSOs makes 
it impossible to definitively test whether these powerful in- 
flows/outflows are found only in the sub-mm detected QSOs 
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Fig. 7. The relationship between Far-Infrared (FIR) luminosity 
and the mass accretion rate onto the black hole for sub-mm de- 
tected QSOs (filled circles with error bars). The equivalent axes 
of star formation rate (SFR) and bolometric luminosity are in- 
dicated on the right and top of the plot, respectively. The sub- 
mm QSOs are close to the highest end of the relationship be- 
tween the two quantities found t o be valid over 6 orders o f mag- 
nitude (black line in the figure; Trakhtenbrot & Netzeri 2010). 
QSOs at z~2 observed in the mid-IR by Spitzer dLutz et aU 
l2008i) are plotted as squares; filled symbols are used for the 
QSOs in common with our sample. 



and not in the sub-mm non-detected QSOs from the parent 
samples. 

Black hole mass estimates used in this work are based on 
the FWHM of the Ha line (FWHMea) and the continuum lu- 
minosity. Using the Ha luminosity instead of continuum lumi- 
nosity lowers the masses by ~3. Using the width of the very 
broad Gaussian component (i.e. the VBLR only instead of both 
VBLR and ILR) to estimate the black hole mass increases 
the black hole mass estimates by ~5 and changes the inter- 
pretation of the black hole growth phase in the QSO sample 
(Fig. |6]l. This uncertainty, coupled with the uncertainty of ap- 
plying reverberation-mapping-based scaling relations at these 
high continuum luminosities limits the reliability of our results, 
and those of other similar works. 

We find no significant differences in the distributions of 
black hole masses or emission-line velocity offsets between the 
subsamples of sub-m m detected and non-de tecte d QSOs from 



the pa rent samples of Pri ddev et al 
(2003). Furthermore, Pridd ev et al 



2003ah and 
l2003al) and 



Omont et al. 



Omont et al 



(l2003h did not find any other predictor of sub-mm loudness. 
The lack of rest-frame optical spectroscopy of the sub-mm non- 
detected QSOs limits any comparison, though it is reasonable 
to assume that their properties are siniilar to othe r z~2 QSOs 



(e.g. t he samples of IShemmer et al.l (120041) and iNetzer et al. 
(l2007h '). 



The sub-mm detected QSOs studied here have the high- 
est rest-frame optical continuum luminosities among sub-mm 
detected galaxies thus far studied in rest-frame optical spec- 
troscopy (e.g. Fig. |5j. These sub-mm detected QSOs host 
black holes in the typical range of luminous QSOs at red- 
shift z~2, and are accreting at close to their Eddington limit 
(Fig. |6ll. In black hole mass and Eddington fraction they are 
indistinguishab l e frorn the optically-luminous QSO sample of 



Shemmer et al. 



(2004). The latter sample has not been stud- 
ied in the sub-mm though it is reasonable to assume that it is 
dominated by sub-mm dim QSOs. If black hole masses in the 
sub-mm detected QSOs are estimated from the VBLR width 
only, then black hole masses increase, and Eddington fractions 
decrease, by a factor of ~5. In this case the sub-mm detected 
QSOs host relatively massive black holes which are accreting 
at a fraction of the Eddington rate. In this case we are un- 
ab le to make a m e aningf ul co mparison to t he Q SO samples 



of Shemmer et al 



(12004 and iNetzer et al.1 (EoO?) as VBLR 
widths are not quoted in their studies. In brief, the sub-mm de- 
tected QSOs are, thus far, indistinguishable from the typical 
z~2 QSO population apart from their relatively high sub-mm 
fluxes. 

INetzer et al.l tOOlh and 



Trakhtenbrot & Netzeri (120101) 



have demonstrated a correlation between SFR and mass accre- 
tion rate in active galaxies, valid over six orders of magnitude. 
The extreme SFRs and mass accretion rates of the sub-mm de- 
tected QSOs places them in the vicinity of the extreme end of 
this trend (Fig. |7]i. However, they are off'set from the relation- 
ship in having a factor ~ 10 higher SFR. Within the sub-mm de- 
tected QSO sample, there is a weak indication that the strongest 
star forming QSOs have the lowest black hole masses and mass 
accretion rates (Fig. |7]l, but this weak trend is w ithin the scatter 
seen when other QSOs are in cluded in the plot. ^ Maiolino et al 



([2007) and lLutzetal] (|2008|) have studied the SFR accretion- 
rate correlation in a sample of optically-luminous QSOs ob- 
served with Spitzer and noted that the SFR potentially flattens 
at the highest bolometric luminosities. Adding the data from 
our sample of QSOs (which changes some of the previous large 
upper-limits on the SFR into known lower values) appears to 
support their findings. 
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Table 2 

The Submm z~2 QSO Sample: Results of the Near-IR Spectroscopy 



Name 


R.A. 


Dec. 


z 


z 


Aobs (A) 


Aobs (A) 


Flux 


Flux 


FWHM 


FWHM 


FWHM 


FWHM 


Flux 




(J2000) 


(J2000) 


rest UV 


Ha 


Ha 


Ha 


Ha 


Ha 


intermed. 


broad 


Lorentz 


Hq 


6000A 












intenned. 


broad 


intermed. 


broad 


(kms"') 


(kms"') 


(kms"') 


(kms"') 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


LBQS B0018-0220* 


00 21 27 


-02 03 34 


2.596 


2.617 


23745 


23745 


5.9E-15 


9.1E-15 


2000 


5000 


2800 


2646 


1.18E-16 


HS B0035+4405* 


00 37 52 


+44 21 33 


2.71 


2.738 


24540 


24500 


l.lE-14 


5.9E-14 


1500 


7500 


3000 


2882 


2.04E-16 


[HB89I 0933+733 


09 37 48 


+73 01 58 


2.528 


2.536 


23210 


23210 


4.0E-14 


5.7E-14 


2700 


6800 


3600 


3504 


2.39E-16 


HSB 1049+4033" 


10 51 58 


+40 17 36 


2.171 


2.189 


20935 


20935 


1.9E-14 


3.4E-14 


2500 


7500 


3500 


3394 


2.20E-16 


HSB 1103+6416 


11 06 10 


+64 00 08 


2.201 


2.216 


21115 


21115 


2.7E-14 


6.6E-14 


3200 


8500 


4500 


4740 


3.39E-I6 


HSB1141+4201* 


1143 52 


+41 45 20 


2.12 


2.152 


20690 


20690 


2.1E-14 


4.6E-14 


2700 


8500 


3900 


3856 


2.62E-16 


SBSB1408+567 


14 09 55 


+56 28 27 


2.562 


2.579 


23493 


23483 


4.6E-15 


1.4E-14 


1800 


4500 


2800 


2854 


8.84E-17 


SBSB 1542+541 


15 43 59 


+53 59 03 


2.371 


2.381 


22193 


22153 


2.3E-14 


5.2E-14 


5400 


9000 


6500 


7120 


1.85E-16 


HSB161 1+4719 


16 12 39 


+47 1 1 58 


2.35 


2.393 


22273 


22273 


7.4E-15 


2.0E-14 


2000 


6000 


3000 


3036 


1.06E-16 


J164914.9+530316 


16 49 14 


+53 03 16 


2.26 


2.274 


21493 


21493 


2.3E-14 


3.1E-14 


2200 


7000 


2800 


2790 


1.75E-16 



Note. — Columns are: (1) QSO name; (2) and (3) position in J2000 coordinates; (4) redshift derived from rest-frame UV lines (references listed below); (5) redshift derived from 
the line center of the intermediate width (ILR) Gaussian component of the double Gaussian fit to the Ha profile. The typical error is about 0.0015 (see below); (6) and (7) observed 
wavelengths of the line centers of the intermediate width (ILR) and very broad (VBLR) Gaussian components fitted to the Ha line, respectively; typical errors are lOA and 20A for the 
inlermediale width and very broad Gaussian components, respectively; (8) and (9) fluxes of the intermediate width and very broad Gaussian components filled to the Ha line, respectively; 
(10) and (11) FWHM of the intermediate width and very broad Gaussian components fitted to the Ha line, respectively; (12) FWHM of the single Lorentzian component fitted to the Ha 
line; (13) FWHM of the Ha line as measured from the summed (ILR+VBLR) double Gaussian fit to the Ha line, and used for the black hole mass estimation; (14) continuum flux near 
6000A, with a K-correction of (1+z) applied. 



Note. — Wavelengths reported are vacuum wavelengths. The continuum fluxes are in units of erg s"' cm"^ A"', and have been K-correcled by factor (1+z). Ha line fluxes are in units 
of erg s"^ cm"-. Errors in fluxes are ^^10- 20% (see text). An asterisk after the QSO name indicates a 850 fim detection (Priddey el al. 2003b). A double asterisk after the QSO name 
indicates that the object is delected at 1.2 mm (Omont et al. 2003) but not detected at 850 fim (Priddey et al. 2003b). The first two QSOs in the table were observed in September 2002 
with the HK grism and the other QSOs in March 2002 with tiie Kb grism. Rest-frame UV ledshifts are from Hewett et al. (1995; LBQS B0018-0220), Hewitt & Burbidge (1989; [HB89] 
0933+733), Siepanian et al. (2001; SBS B1408+567), SDSS (SBS B1542+541, HS B1049+4033 and HS B1103+6416), and from flie Hamburg Bright Quasar Survey (remaining four 
QSOs; Hagen et al. 1999; Engels et al. 1998). 
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Table 3 

The Sub-mm z~2 QSO Sample: Luminosities, Black Hole Masses, and Eddington Rates 



Name 


log i-M 


log Lfir 


logMfi//(//a) 


logMfJH(L5100) 


logL/LEM 




(ergs s"') 


(i.0) 


(Mq) 


(Me) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


LBQS B0018-0220 


47.30 ± 0.04 


13.17 ±0.76 


8.62 ±0.14 


9.07 ±0.10 


0.06 ± 0.06 


HS B0035+4405 


47.58 ± 0.04 


12.90 ± 0.84 


9.17 ±0.14 


9.29 ± 0.09 


0.08 ± 0.10 


[HB89J 0933+733 


47.57 ± 0.04 


12.93 ± 0.80 


9.29 ±0.14 


9.45 ± 0.09 


-0.10 ±0.09 


HS B1049+4033 


47.39 ± 0.04 


12.88 ± 0.79 


9.06 ±0.14 


9.33 ± 0.09 


-0.14 ±0.07 


HS 81103+6416 


47.59 ± 0.04 


12.96 ± 0.83 


9.52 ±0.14 


9.73 ± 0.08 


-0.36 ± 0.08 


HS B 1141+4201 


47.44 ± 0.04 


12.88 ± 0.84 


9.23 ±0.14 


9.48 ± 0.09 


-0.23 ± 0.08 


SBS B 1408+567 


47.16 ± 0.04 


13.37 ± 0.65 


8.79 ±0.14 


9.06 ± 0.09 


-0.06 ± 0.08 


SBSB 1542+541 


47.40 ± 0.04 


12.94 ± 0.84 


9.87 ±0.14 


10.00 ± 0.08 


-0.79 ± 0.09 


HSB161 1+4719 


47.16 ±0.04 


13.02 ± 0.74 


8.87 ±0.14 


9.12 ± 0.09 


-0.12 ±0.08 


J164914.9+530316 


47.32 ± 0.04 


13.03 ± 0.76 


8.88 ±0.14 


9.13 ± 0.09 


0.02 ± 0.08 



Note. — Columns are: (1) QSO name; (2) Log of the bolometric luminosity and its la error; (3) Log of 
the FIR luminosity (calculated from the sub-mm flux following McMahon et al. (1999) with = 1.5 and Tj 
= 40 K) and its Icr enm (which reflects only on the error of the submm flux). (4) Log of the black hole mass 
as estimated from LHa and FWHMna, and its la error (see text); (5) Log of the black hole mass as estimated 
fromLsioo and FWHMna, and its Icr error (see text); (6) Log of the accretion rate as a fraction of the Eddington 
accretion Rate. 



Note. — The flux to luminosity conversion was made using a luminosity distance derived within the standard 
cosmology with parameters listed in Sect. 1. The errors listed in this table reflect measurement errors and some 
of the declared errors in the scaling relationships used. More realistically, the scaling relations used to estimate 
black hole masses are expected to have errors of a factor of a few. 



